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Abstract: We investigate the effect of a magnetic field on the vacuum of low- 
temperature QCD at large-iV c in the presence of a chiral chemical potential, using 
the holographic Sakai-Sugimoto model. Above some critical chemical potential we 
find an instability, which triggers a decay of the homogeneous vacuum to a non- 
homogeneous configuration with a spiral form, which we construct explicitly. We find 
that this decay is suppressed for sufficiently large magnetic field, and determine the 
critical strength. We find that the new vacuum does not exhibit the chiral magnetic 
effect. This is related to the way the chiral chemical potential is introduced. We 
discuss an alternative way of introducing the chiral chemical potential that leads to 
a nonzero chiral magnetic effect. 
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1 Introduction and discussion 

Recent work by various authors has provided evidence that in QCD, a magnetic field 
leads to several interesting phenomena, for values of the magnetic field which are 
substantially smaller than anticipated in earlier studies. The origin of these effects 
lies in the chiral axial anomaly and the additional couplings it provides. For instance, 
it has been shown [1] that the axial anomaly is responsible for the appearance of 
pion domain walls, which carry baryon charge. These are stable for sufficiently large 
baryon chemical potentials. Importantly, they form for relatively small critical values 
of the magnetic field, as this critical value scales with the pion mass. 
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In the presence of an axial chemical potential [2, 3], an external magnetic field 
can trigger the appearance of a vectorial current in the direction of this field, 



(1.1) 



The origin of this so called chiral magnetic effect (CME) is the asymmetry between 
left- and right-handed fermions (parametrised by the axial potential (ia), which leads 
to separation of electric charge. Related to this is the chiral separation effect [4-6], 
where instead of an axial potential, a baryonic quark potential is introduced. The 
magnetic field now leads to a separation of chiral charge, and an axial current is 
generated, 



Combining these two effects leads to evidence for the existence of a so-called chiral 
magnetic wave [7]. All these effects suggest that the consequence of a magnetic field 
may be much more important (and potentially easier to see in experiments) than 
previously thought. 

These effects were originally derived using a quasi-particle picture of chiral charge 
carriers, which is only valid at weak coupling. They have meanwhile also been 
confirmed in a number of holographic descriptions of the strong coupling regime. The 
pion gradient walls were found in the Sakai-Sugimoto model in [8, 9] (see also [10]), 
and there are even abelian analogues of this solution [11] which involve an rj- meson 
gradient (and of course no baryon number in this case). Whether or not holographic 
models exhibit the CME or CSE is subject of more controversy, as it requires a 
careful definition of the holographic currents. An analysis of the D3/D7 system was 
given in [12, 13] and we will comment extensively on the Sakai-Sugimoto model later 
in this paper. 

A major question is what happens when the chemical potentials are large enough 
so that they trigger a condensation of bound states. Even without a magnetic field 
and without an axial anomaly, various models predict the formation of a non-isotropic 
(although homogeneous) vector meson condensate, once the axial chemical potential 
becomes of the order of the meson mass (see [14] for an analysis in the Sakai-Sugimoto 
model relevant here, and a list of references to other works). When the axial anomaly 
is present, the condensate is typically no longer even homogeneous, but forms a spiral 
structure [15]. A lot of work in this direction focuses on high temperature models, 
but in fact such condensation already happens at low and zero temperatures [14, 15] . 

In the present paper, we therefore examine the Sakai-Sugimoto model in the 
presence of a magnetic field and an abelian chiral chemical potential, and set out to 
determine whether there is an instability against decay to a non-homogeneous chiral 
spiral also in this case. We will also analyse whether the ground state exhibits a 
chiral magnetic effect, and whether there is an r/-gradient in this case. We will only 
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consider the confined chirally broken phase of this model, i.e. the low-temperature 
behaviour. 

Our main results are as follows. First, we determine the location of the phase 
boundary between the homogeneous solution of [8, 9] and a new chiral spiral like 
solution 1 . We construct the non-homogeneous chiral spiral like solution explicitly. 2 
We establish that the latter still does not exhibit the CME, but point out that 
this conclusion relies crucially on subtleties related to the precise way in which the 
chemical potential is introduced. Finally, we show that a magnetic field tends to 
stabilise the homogeneous solution, and there is a critical value beyond which the 
chiral spiral ceases to exist. The parameter space in which we need to scan for 
non-linear spiral solutions is rather large and the numerical analysis is consequently 
sometimes tricky; we comment on details of this procedure in the appendix. 

An important issue in this analysis is the precise way in which the chemical 
potential is included in the holographic setup. It was recently pointed out [13, 16] that 
the usual approach, which roughly amounts to identifying the chemical potential with 
the asymptotic value of the gauge field, may be incorrect in case the associated charge 
is not conserved (like, in our case, the axial charge). However, we can define \xa as the 
integrated radial electric flux between the two boundaries of the brane. Following 
the field theory arguments of [16] we will argue that there are then two natural 
formalisms for introducing the chemical potential. In formalism A we introduce lia 
through boundary conditions in the component A§ of the gauge field. In formalism 
B the chemical potential instead sits in A z . In the presence of the axial anomaly 
these two formalisms are inequivalent. 

In formalism A, we find that the non-homogeneous phase characterised by wave 
number k, is dominant and there is a preferred value of this wave number k = fc g . s .. 
When the magnetic field is small compared to the chiral chemical potential B <C fx a, 
kg. B . depends only very weakly on the value of lla- This result is consistent with 
our previous work [15]. For sufficiently large magnetic field we find that the non- 
homogeneous ground state is suppressed. In formalism B, on the other hand, the 
homogeneous state always has lower energy than the non-homogeneous one, and 
there is hence no chiral magnetic spiral. 

As far as the chiral magnetic effect is concerned, we find that in formalism A 
the effect is absent. This result is consistent with previous calculations on the Sakai- 
Sugimoto model [17] . In formalism B, on the other hand, there exists a non-zero chiral 
magnetic effect. This seems to be more in line with recent lattice calculations [18], 

■"■The stability analysis of the homogeneous solution presented in [11] is incorrect in some crucial 
aspects; we will comment on the details and the consequences later. 

2 The chiral magnetic wave mentioned earlier is different from the chiral magnetic spiral, because 
the chiral magnetic wave involves the density and the component of the current parallel to the 
magnetic field, while the chiral magnetic spiral involves the components of the current transverse 
to the magnetic field. 
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chiral model calculations [19] and holographic bottom-up models [20, 21] for the 
confined chirally broken phase of QCD. However, a main shortcoming of formalism B 
is that the inhomogeneous phase, whose existence is indicated by our perturbative 
analysis, is absent in the full nonlinear theory. Since the perturbative analysis is blind 
to the subtle issues on how one introduces the chemical potential in the holographic 
setup, we are inclined to think that, as it stands, formalism B is incorrect. Whether 
formalism A is correct on the other hand, or also needs to be altered, is an open 
question at the moment, and needs a separate study, which we leave for future work. 



2 Chemical potentials, currents and anomalies 

2.1 Effective five-dimensional action 

In order to set the scene and to introduce our conventions, let us start by giving a 
brief review of the basics of the Sakai-Sugimoto model [22, 23]. For a more detailed 
description of the features of this model which are relevant here, we refer the reader 
to the original papers or to [15]. 

The Sakai-Sugimoto model at low temperature consists of Nf flavour D8-branes 
and Nf anti-D8-branes which are placed in the gravitational background sourced 
by a large number N c of D4-branes which are compactified on a circle of radius R. 
In the simplest set up, the probes are asymptotically positioned at the antipodal 
points of a circle, while in the interior of the geometry they merge together into 
one U-shaped object. The gauge theory on the world-volume of the probe brane 
is nine-dimensional, containing a four-dimensional sphere, the holographic direction 
and four directions parallel to the boundary. By focusing on the sector of the theory 
that does not include excitations over the S A , one can integrate the probe brane 
action over this sub-manifold and end up with an effective five dimensional DBI 
action on the probe brane world-volume. By expanding this action to the leading 
order in the string tension, one ends up with a five dimensional Yang-Mills theory 
with a Chern-Simons term. For Nf = 1 we can write the action as 

S = S YM + S CS = -^J d*xdz^F mn T mn + j e lmn ™ J d 4 xdzA e F mn F pq . (2.1) 

where the indices are now raised or lowered using the effective five dimensional metric 
g mn . This metric is defined as 

d4) = g mn dx m dx n = M^ K K 2 J 3 T^d^dz" + K'^dz 2 , (2.2) 

where K z = (1 + z 2 ). The x M directions are parallel to the boundary and z is the 
holographic direction orthogonal to the boundary. The coupling constants k and a 
are given in terms of the number of colours N c , the compactification mass scale Mkk 
and the 't Hooft coupling A by 

_ Vrtg s N 2 M KK _AA^_ _ _A^_ 

K 108tt 2 216^' a 24tt 2 ' [ ] 
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The Chern-Simons term written in (2.1) is valid only for a single D8 probe, i.e. for 
a U(l) gauge theory on the brane world- volume. 



2.2 Symmetries and chemical potentials 

Holographic models encode global symmetries of the dual gauge theory in the form of 
gauge symmetries in the bulk theory, and these relations hold both for the closed as 
well as the open sectors of the string theory. The same is true for the Sakai-Sugimoto 
model at hand, where there are two independent gauge theories living near the two 
boundaries of the flavour D8-D8 brane system. These U(Nf)i, and U(Nf)n gauge 
symmetries correspond to the global U(Nf) x U(Nf)n flavour symmetries of the dual 
gauge theory. In the low-temperature phase of the Sakai-Sugimoto model that we 
are interested in, the two branes are connected in the interior of the bulk space, and 
thus the gauge fields Am and A M are limits of a single gauge field living on the two 
connected branes. Therefore, one cannot independently perform gauge transforma- 
tions on these two gauge fields, but is constrained to gauge transformations which, 
as z — > ±00, act in a related way. Specifically, since near the boundary a large bulk 
gauge transformation acts as Al/r — > gL/RAL/Rg~[j R , then clearly any state (and in 
particular the trivial vacuum A = 0) is invariant under the vectorial transforma- 
tions, i.e. those transformations for which gi = gR- This means that the vector-like 
symmetry is unbroken in this model. On the other hand, the fact that the branes 
are joined into one U-shaped object means that the axial symmetry is broken, which 
corresponds to the spontaneous breaking of the axial symmetry in the dual gauge 
theory. The corresponding Goldstone bosons can be seen explicitly in the spectrum 
of the fluctuations on the brane world-volume. 

The relationship between the bulk gauge field and the source and global symme- 
try current of the dual gauge theory is encoded in the asymptotic behaviour of the 
former. More precisely, the bulk gauge field Am(z,x^) behaves, near the boundary 
and in the A z = gauge, as 



Here p u parametrises the normalisable mode, while a !/ (x M ) describes the non-norm- 
alisable behaviour of the field. The latter is interpreted as a source in the dual field 
theory action, where it appears as 



Hence the expectation value of the current corresponding to the global symmetry in 
the gauge theory is given by 



A v (x», z) ->• a v (x») (l + 0(z- 2 ^ + p v {xr)-(l + 0(^- 2/3 )) 



(2.4) 




(2.5) 



J"{x)± 



5S 



(2.6) 



SA^(x, z —7- ±00) 
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When the bulk action is just the ordinary Yang-Mills action (in curved space), this 
expectation value is the same as the coefficient p v in the expansion (2.4). However, 
in the presence of the Chern-Simons term, the coefficient p v is different from the 
current, as we will explicitly demonstrate for the system at hand in section 2.3. This 
difference has been a source of some of confusion in the literature. Its importance 
has recently been emphasised in the context of the chiral magnetic effect in the 
Sakai-Sugimoto model in [24]. 

From (2.5) we also see that adding a chemical potential to the field theory cor- 
responds to adding a source for J°, which implies the boundary condition for the 
holographic gauge field A v {x) = p5„o- 

For the Sakai-Sugimoto model, the bulk field A m living on the D8-branes has two 
asymptotic regions, corresponding to each brane, and hence there are two indepen- 
dent chemical potentials //£ and fiR which can be separately turned on. Instead of 
left and right chemical potentials one often introduces vectorial and axial potentials, 
defined respectively as \ib = \{^r + and fx a = \{hr — Ml)- For the Nf = 2 
case, the vectorial and axial chemical potentials for the U(l) subgroup of the U(2) 
gauge group on the two D8-branes correspond to the baryonic and axial chemical 
potential in the dual gauge theory, while the non-abelian £77(2) chemical potentials 
are mapped to the vectorial and axial isospin potentials. 

In what follows we will also be interested in studying the system in the presence 
of an external (non-dynamical) magnetic source, which will be introduced by turning 
on a nontrivial profile for the non-normalisable component a v {x) of the bulk field. 

2.3 The Chern-Simons term, anomalies and the Bardeen counter term 

The symmetries and currents of the model discussed in the previous section are, 
however, valid only at leading order in A -1 . At the next order in A -1 , the Yang-Mills 
action on the brane world- volume receives many corrections, among which the Chern- 
Simons term. This term turns out to be crucial for the existence of a nontrivial ground 
state of the system in the presence of the external magnetic field, a situation which 
we will study in the following sections. On manifolds with boundaries, however, the 
Chern-Simons term is manifestly gauge non-invariant and it also spoils conservation 
of the vectorial and axial currents. Both of these reflect the fact that the Chern- 
Simons term is the holographic manifestation of the vector and axial gauge anomalies 
in the dual theory. 

Apart from the Chern-Simons term there are of course also various other correc- 
tions. One important class of terms is coming from the expansion of the DBI action. 
In previous work [15] we have seen that the qualitative picture of chiral spiral for- 
mation at vanishing magnetic field does not change when these corrections are taken 
into account. We will assume that this holds true here as well. For the other higher 
derivative corrections, it is important to note that we will find that our ground state 
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has a momentum scale of order Mkk> not l/l s , so higher derivatives can typically be 
ignored as long as M KK l s is small. 

Returning to the role of the Chern-Simons term in describing anomalies, let us 
start by decomposing the gauge field in terms of the axial and vectorial components, 

as 

•Afn^X) z) A m {x 1 z) -\- A m {x , z} . (2-7) 

These two components transform under under an inversion of the holographic coor- 
dinate z — > — z as 

A v J A {-z,x) = ±A v J A (z,x) , A v J\-z,x) = TA v J A {z,x) . (2.8) 
Furthermore, the /i component of these fields are related to the dual gauge fields as 
Al(z = ±00, x) = Al(x) , A A (z = ±00, x) = T A A (x), (2.9) 

where the non-calligraphic A^ A (x) are the boundary vector and axial- vector gauge 
fields. When written in terms of these vectorial and axial potentials the action reads 



S = J d' X( iz{ - K -^g \iy n Tl n + n n ?L\ 

+ pq J~ mn J- pq + Ai (J~ mn J- pq + J~ mn J~ pq )^ I . (2.10) 

We now want to compute the currents following the prescription (2.6). The variation 
of the action can be written as 



5S= d 4 xdz 



dAY m v 



5AJ + 



Ui ~ _ a sprat 

dA A rn 



5Ai 



+ d m [vfdAj + Vf6A A ] } , (2.11 



where the derivatives of the Lagrangian are given by 



dC 
dAj 

dC 
dAf 

'rym£ 
' V 



V 



ml 



a 



Imnpq ~pA -pV 
ran'' pq J 



^ Imnpq ( -r-V pV 



\" mn^ pq 1 " ran" pqJ ' 



pq) 



dC 



d(d m AY) 

dC 

d(d m A A ) 



-2kV=W + ae^{AlF A + A A T v vq ) , 
-2K^r g Tf + ae m ^\A v n T v vq + A A n T A q ) . 



(2.12) 
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Note that V$f A are antisymmetric in m -H- £. Imposing that the bulk term in the 
variation (the first line in (2.11)) vanishes gives the equations of motion, 



(2.13) 



We will need these shortly to show that the currents are not conserved. The boundary 
term in the action variation can be written as 



5S } 



bdy 



d 4 xd^ \d z [P z /8Al + V Z /SA$] + d. 



PifSAY + PfSAi 



(2.14) 



This implies that the holographic vector and axial currents are given by the expres- 
sions 



J v (x) = -Ak lim y=gjy] - 2ae^°{A v v F* + A$F f 
J%x) = Ak lim [V=0J7] - 2ae^{A v v Fl + A A v F t 



pa) i 



At?A\ 
pa) i 



(2.15) 



where we have used the boundary conditions (2.9) and F P V J A = d^AX — d v A V J A are 
the boundary field strengths. 

Using the equations of motion (2.13) we can explicitly show that these two 
currents are not conserved due to the presence of the Chern-Simons term. One finds 



d^Jyix) = 2 lim «9 M [^gVf \ = -2 lim 



-9 



dC 
dAj 



" c fiu pa ) 



dMx) 



-2 lim \\l ZI gV z /\ = 2 lim 



- dC 



(2.16) 



a 
2 



pupa/rpV rpV _,_ t?A t?A\ 
\ p.v pa 1 pv pa) i 



We see two things here. Firstly, the anomaly (i.e. the right-hand side of above 
equations) is indeed sourced by the Chern-Simons term. Secondly, the anomaly is 
present only if the boundary value of the bulk gauge field strength is non-vanishing. 
In other words, only in the presence of an external field does the anomaly show up. 

While the axial anomaly is not problematic (it just reflects the fact that this 
symmetry is no longer present in the dual quantum field theory), this is not the case 
with the vectorial symmetry. In QED coupled to chiral fermions one has to require 
that the vector current is strictly conserved, since its non-conservation would imply 
a gauge anomaly. It is possible to make the vectorial current conserved by adding 
extra boundary terms to the action, as was first shown by Bardeen in [25]. 

In the holographic setup, one expects that such a Bardeen-type counter term will 
appear, but this type of term should come from the requirement that the full theory 
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in the bulk in the presence of the boundary is gauge invariant under vectorial gauge 
transformations. Let us therefore consider a generic vectorial gauge transformation 
in the bulk 

8Aj = d e A v (x, z), 5Af = , (2.17) 

where A-y(x,z) is a function even in z. Under this transformation, the action (2.10) 
is invariant up to a boundary term, 

= J d 4 xdz(d e A v )d m Vf 

= \ t emnm f d 4 xdz(d e A v )^ n ^ g (2.18) 
= ae tmn ™ [ d 4 xdzd m l(d,A v )At^ q ] . 



Therefore, if we want to impose invariance of the action under (2.17) we need to add 
the anomaly counter term correction 

^n. = -ae" mnm j d A xdzd m [AjA^ q ] 

(2.19) 

_ « Jmnpg f j4 j \_aV -pA -pV , a A -pV pVl 

We see that there are two contributions of this surface term: at the "holographic" 
boundaries z — > ±oo and at the boundary at spatial infinity \x\ — > oo. Its contribu- 
tion at holographic infinity is indeed the Bardeen counter term as derived in quantum 
field theory [25], and which was, in the holographic setup, initially postulated (added 
by hand) in [24]. We see that its presence automatically follows from the requirement 
of classical gauge invariance of the bulk theory in the presence of the boundary. 

The contribution at spatial infinity, would typically vanish, as all physical states 
in the system are localised in the interior. However, in the presence of external sources 
which generate an external magnetic field that fills out the whole four-dimensional 
space, like the one we will be considering, it is not a priori clear if this is true, and 
one has to be careful about possible extra contributions to the action and currents. 
It will however turn out that for our non-homogeneous ansatz, these extra terms are 
irrelevant and that only the Bardeen counter term is non-vanishing. 

Let us continue by showing the effect of adding the S an _ term on the currents. 
The total action now reads 

S = S + S^n. = J dVdz j — —\/—g \Ty n J\ nn + J^^-^mn] 

+ ^e im ^Ai( 3 Tl n T v m + } = j d'xdzC . (2.20) 
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The variation of this new action can again be written in the form (2.11), with C and 
P, instead of C and P, i.e. 



dC 



0. 



dAj 
dC a 



> 



Imnpq / o -pV -pV . 77 A 77 A •> 



(2.21) 



The equations of motion are of course unchanged (as the Bardeen counter term is 
only a surface term), while the new currents are obtained as 

bAfj, 1 (x) oAp (x, z = oo) z ^°° 
Explicitly, the expressions read 

fy(x) = -4k lim [V=gJ^] ~ Qae^A^F^ , 

z— >oo ' 

Jftx) = Ak lim [y/=j7*f] - 2ae^A?F* . 
The divergences of these currents are 



(2.23) 



d (i 4(x) = o, 



(2.24) 



where we used again the equations of motion (2.13). This clearly shows that the 
vector current is now conserved, while the anomaly is seen only in the axial sector. 
When the coefficient a is taken from string theory, the non-conservation of the axial 
current is exactly the same (including the numerical factor) as in QED coupled to 
external fermions [26] . In what follows we will work with these renormalised currents 
and action. We should also emphasise that when one considers chemical potential 
for a symmetry which is anomalous (as it is case here), just knowing the corrected 
action (2.20) may not be enough. Also, one has to be careful about the boundary 
conditions one has to impose on the states, as not all boundary conditions are allowed. 
We discuss these subtle issues in section 3.4. 



-10- 



2.4 The corrected Hamiltonian 



The Lagrangian density of the action after inclusion of the anomaly counter term 
can be written as 



C 



-/C-\ 



pOa pV I pab pV , x-Oa pA , 1 -r-ab -r-A 
V Oa ' V r - r ab~ t ~- r A- r 0a~r r,-' A - r ab 



a 
+ 4 £ 



Oabcd 



^(3^K + ^^) + 4^(3J- ^ + J-^) . (2.25) 



The conjugate momenta associated with the vector and axial gauge fields thus take 
the form 

^Oa q,„ / ~ -pOa I q „ r 0abcd a A pV 



d(doAV) 
A d(d A£) 



-2m 



(2.26) 



P 



Oa 



We then obtain the on-shell Hamiltonian as H = i^Buik + H^ y , where the two 
contributions read 



Bulk 



K j d 3 xdzy/—g 
k l d 3 xdzy/—g 



pOa pV I pab pV pOa pA , ^ pab pA 
V Oa ' 2^ V Alb ^ A ^Oa "T 2^ A ''ab 



1 



X"Ua X" 1 pab p 
—J~ •> Oa + ^ a6 



(2.27) 



Bdy 



d 3 xd^ d„ 



W v Al + fl a A A, 



a-^o 



(2.28) 



Here we have used the gauge field equations (2.13) for the time component 
(generalised Gauss law). 



3 The spatially modulated phase 

Having settled the issue of how to deal with the Chern-Simons term in the presence of 
external fields for the Sakai-Sugimoto model, we now want to find the ground state 
of the system, at strong coupling, in the presence of external magnetic field and 
non-vanishing axial and baryon chemical potentials. Based on the weak coupling, 
partonic, arguments which were mentioned in the introduction, we expect that the 
ground state should be a chiral magnetic spiral like configuration. In particular, we 
expect it to be non-homogeneous. 

So far, for Sakai-Sugimoto model a non-homogeneous ground state was con- 
structed in the presence of large enough axial chemical potential but with no exter- 
nal fields in [15] (at low temperature; see [27] for a high-temperature analysis). The 
main reason why such a state appeared was due to the nontrivial Chern-Simons term. 
We now want to see if this state persists and how it is modified once the external 
magnetic field is introduced. 
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3.1 Magnetic ground state ansatz in the presence of chemical potentials 

We are interested in studying the ground state of the system at non-zero axial chemi- 
cal potential Ha and non-zero baryon chemical potential Hb- 111 addition, we will turn 
on a constant magnetic field in the x l direction B = Bx±. The boundary conditions 
associated with this physical scenario are 

A(z = ±00) = -B x x, Aq(z = ±00) = THa + H-b ■ (3.1) 

We will see that in fact our ansatz is insensitive to the baryon chemical potential, 
but we will keep it in the formulas for a little while longer. 

Let us now consider our particular ansatz, in the gauge A z = 0. First, we want 
to introduce the vectorial and axial chemical potentials, hence we turn on Ao = f(z) 
with the above boundary conditions. Second, we need to introduce the (constant) 
magnetic field in the boundary, hence we turn on a component of A transverse to the 
direction of B: A] 3 (x 2 , x 3 , z). In principle this function could depend on z, but the 
Bianchi identity tells us that in fact the magnetic field is independent of z. Therefore, 
we have for A^ 

1 -> B 

A B (x 2 , x 3 ) = -B x x = — [-X3X2 + x 2 x 3 ] . (3.2) 

Next by looking at the equations of motion for the A B component parallel to the 
magnetic field, we see that this component also has to be turned on, and is only a 
function of z, 

A L B = a(z) Xl , (3.3) 

Finally, we expect that a chiral magnetic spiral will appear in the direction transverse 
to the external magnetic field 

A w (xi,z) = h{z) [cos(A;a;i)a;2 — s\n(kxi)x 3 \ , (3.4) 

i.e. it represents the chiral wave transverse to the boundary magnetic field satisfying 

VxA^ v = kA^ V) (3.5) 

where k = ±\k\ and k is the spatial momentum. So in summary our ansatz is given 
by 

J% = fv(z) , A} = f A (z), 
-> B 

Ay = — [-X3X2 + a; 2 x 3 ] + h v (z) [cos(kxi)x 2 - sin(/cxi)x 3 ] + a v (z)x 1 , (3.6) 
Aa = cla{z)xi + h A (z) [cos(kxi)x 2 — sm(kxi)x 3 ] , 
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where the fields satisfy the boundary conditions 

f V (z ->• ±00) = Hb , }a{z ± 00) = =F/iA ■ 

hy(z — > ±00) = , h,A(z — > ±00) = , 

ay{z — > ±00) = , o,a{z — > ±00) = =Fj . 



(3.7) 



The Gauss law, i.e. the zeroth component of the equation of motion (2.13), 
is automatically satisfied for our ansatz. The remaining equations reduce, after 
integrating one of them with integration constant p, to 



V-gg zz g m d z f = 3 



a 



Ba + -h 2 
2 



d z [V=gg zz g xx d z a] + 3-BdJ = , 

d z [V=jj g zz g xx (d z h)} - ^Tg {g xx fk 2 h + 3-kdJh = . 

K 



(3.8) 

(3.9) 
(3.10) 



Restricting now to the metric (2.2), and substituting eq. (3.8) into (3.9) and (3.10) 
we obtain our master equations, 



K z d z f 
K z d z 

K z d z 



-b- -kh 2 
2 



K z d z b 
K z d z h 



-B 2 


b+-kh 2 




to 



- K 2/3 k 2 h - kh 



b+-kh 2 



where 



(3.11) 
(3.12) 

(3.13) 
(3.14) 



b = Ba — p , 

and we have also introduced a set of dimensionless variables /, h, a, k, p and B 
defined by 



f = XM KK f , h = XM KK h , a = XM KK a. 
k = M KK k , ~p = \Ml K p , B = \Ml K B, 



(3.15) 
(3.16) 



with A = A/(277r). These coupled equations are in general not solvable analytically, 
except in the special case k = 0, which we will review next. 

3.2 Review of the homogeneous solution 

Before embarking on the full task of finding the non-homogeneous solutions to the 
equations of motion, we will in this section first review the homogeneous solution 
(i.e. the solution for which k = 0) in the presence of a constant magnetic field. This 
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is an abelian version of the solution first constructed in [8, 9], see also [17]. In the 
homogeneous case there is no transverse spiral, i.e. h = and the equations of motion 
simplify to 

dgf =p- Ba, 

/ *■ \ ■ ( 3 - 17 ) 

(dl - 5 2 ) a = -Bp , 

with z = arctan z. These equations can be integrated exactly, and the solution takes 
the form 

a(z) = cosh(I? arctan z) H — ^- sinh(5 arctan z) + 4- , 

B , B , B (3.18) 

C C 

f(z) = — sinh(i? arctan z) — cosh(5 arctan z) + fo , 

B B 

where Ca,Cb,P and fo are four integration constants for the two second order dif- 
ferential equations. The corresponding field strengths take the form 



1 + z 2 

T _ 1 



Ca sinh(5 arctan z) + Cb cosh(5 arctan z 
Ca cosh(S arctan z) + Cb sinh(5 arctan z 



(3.19) 



and Ca/b = ^MkkCa/b- In terms of the baryonic and axial chemical potentials the 
boundary condition on f(z) reads 

f(z -> ±oo) = p L/R = p B =F Pa ■ (3.20) 

From the expression for f(z) these potentials are related to Ca and Cb by 

B B 

Ca= 7 — 7\PA, C B = t — tt-(A-b — fo) , (3.21) 

sinhffm cosh Ub) 

where Pa/b = ^MkkPa/b- In analogy with the analysis of [1] we are here interested 
in configurations in which there is a non-vanishing pion (or rather, r]') gradient in 
the direction of the external field. Since we are working in the A z = gauge, a pion 
field will appear as part of the axial, non-normalisable component of all A^s (see 
[22, 23]). Hence, we impose the boundary condition 3 

a{z ->■ ±oo) = Tj ■ (3.22) 

We should note here that this is not the most general boundary condition for the field 
a(z), since we have set the even part to zero. Having thus reduced the parameter 



3 We stick to this notation, introduced in [17], but want to emphasise that even though using the 
symbol j suggests that the asymptotic value of a is a current, it is not. 
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space to the set {ftA, ftB, j}, we have the relations 



C, 



B 



sinh ( f B 



ftA 



B 



sinh ( f B 



p = -B coth ( -B ) fi A 



TT 



j , /o = Ab - coth -B j 



7T 



(3.23) 



where j = XM^xj- Using these relations we can rewrite /(z) and a(z) as 



ffl u 



sinh (.B arctan z) 

ftB ~ ftA y—^ H j 

sinh f B 



cosh(£? arctan z) 



sinh ( \B 



coth -5 



(I- 



ft A 



cosh(5 arctan z) 



sinh ( IB 



coth ( ^5 



(3.24) 



; sinh (S arctan z) 



sinh 



Since the constant does not appear in the Hamiltonian, it is effectively a free 
parameter, which we are free to set to zero. As expected, we find that the baryon 
chemical potential has no effect on this abelian system. In contrast, minimising 
the Hamiltonian will impose a constraint on the axial chemical potential and j, as 
expected for physical systems. 

In summary, we have two physical boundary values, fx a and j, which are im- 
plicitly expressed in terms of the two parameters Ca and Cg. At any given fixed 
values of lla and j, we will want to compare the homogeneous solution given above 
to possible non-homogeneous condensates and determine which of the two has lower 
energy. 

3.3 Currents for the homogeneous and non-homogeneous ansatz 

The work of [17, 24], which studied the homogeneous solution discussed above, re- 
sulted in the interesting conclusion that there is no chiral magnetic effect present in 
the holographic Sakai-Sugimoto model. There are some subtleties with this which 
were pointed out in [16], to which we will return shortly. However, their result also 
leaves the open question as to whether there is a chiral magnetic effect for more 
general solutions to the equations of motion, for instance the non-homogeneous ones 
which we consider here. 

So even before we find the full non-homogeneous solution, an important lesson 
might be learnt from an evaluation of the corrected holographic currents (2.23) for 
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the non-homogeneous ansatz. For our ansatz (3.6), the corrected currents become 
J v = -4k lim [^g zz g 00 d z f v ] - \2aBj 

z— >oo 

= -12a lim [Ba A + kh v h A ] - UaBj = , 

z— >oo 

J v = -4ac lim \^[ ZI gg zz g xx d z ay\ + 12aBfi A 

z— >oo 

= 12aB lim / A + 12aBfi A = , 

z— >oo 

J v = _4 K ii m \y/=gsf z f x d z jC[\ = -4kM^ k lim [K z d z h v ] cos(k Xl ) , 

z— >oo z— >oo 

Jy = -4k lim \yT I gg zz g xx d z A\\ = 4kM£ k lim [KAV] sin(fcn) , ,„ 9 ^ 

J° A = Ak lim [V=gg zz g 00 dJ A ] 



lim 

z— »oo 



k k 
12a [ Bay + -h\ + -h\ ) - Anp 



-4rep, 

J\ = 4k lim [y/=gg"grd z a A ] = 4kM^ k lim [K z d z a A ] , 

z— >co z— >oo 

= 4k lim \^gg zz g xx d z Af\ = 4kM£ k lim [K z d z h A ] cos(k Xl ) , 

z— >oo z— >oo 

J! = 4k lim \\^gg zz g xx d z Af] = -4kM| k lim [if^/i J smCkxx) . 

z— >oo z— >oo 

where we have used the gauge field equations for the components I = and £ = 1. 

From (3.25) we see two important facts. First, the density of particles carrying 
baryonic charge is zero. This confirms once more that there is nothing baryonic in the 
solutions under consideration, in agreement with the fact that the baryon chemical 
potential \xb decouples completely. 

The second observation is that the component of the vector current in the direc- 
tion of the external magnetic field is zero. This is in sharp contrast to what one would 
expect if there was a chiral magnetic effect present. We should, however, emphasise 
that it is possible to define corrected currents which are different from those above, if 
one decides to deal with the anomalous symmetry in a different way, see section 3.4. 
However, this alternative method, although it produces the chiral magnetic effect, 
suffers from other shortcomings, as we will explain. 

The expressions (3.25) simplify further when we restrict to the homogeneous 
ansatz, for which one obtains 



T° — T 1 — T 2 — T 3 — n 

Jy — Jy — Jy — Jy — U , 

J° A = AkM^ k B fi A coth(^S) = 12aBfi A coth(^B) , 
J \ = ~^Ml K B 3 coth(p) = -12aB(ji B - /„) , 

J 2 A = J 3 A = o. 



(3.26) 
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We should emphasise once more that the vector (baryonic) currents vanish only when 
the contribution from the Bardeen term in the action is properly taken into account. 
One could in principle evaluate an "abelianised" version of the baryon number (sec- 
ond Chern class), ~ f F 3z Fi 2 . However, this expression is strictly speaking valid 
only in a nonabelian system, and in this situation one expects that the charge com- 
puted using the corrected conserved current J v should coincide with this topological 
number. 

Our analysis shows that the fact that there is no chiral magnetic effect in the 
Sakai-Sugimoto model is not due to the simplified homogeneous ansatz, and that it is 
also not a consequence of the details of any numerical solution which we will present 
later; rather, the chiral magnetic effect is absent for the entire class of solutions 
captured by the ansatz (3.6). The contribution of the Bardeen terms required to 
make the vector currents conserved is crucial for the absence of the chiral magnetic 
effect. 

We finally see that, in contrast to the homogeneous solution, where no vector 
currents are present at all, the non-homogeneous system exhibits transverse vector 
currents. This shows some resemblance to the chiral magnetic spiral of [28]. 

3.4 Chemical potentials for non-conserved charges 

Given the rather convincing quasi-particle picture of the origin of the chiral magnetic 
effect at weak coupling [2, 3], it is somewhat surprising that it is not present in the 
model at hand at strong coupling. A reason for this discrepancy has been suggested 
in [16], in which it was emphasised that one should be careful in computing the effects 
of a chemical potential in theories for which the associated charge is not conserved. 

The main observation made in [16] is that there are two ways to introduce a 
chemical potential into a thermal quantum system. One is to twist the fermions 
along the thermal circle, i.e. to impose 



instead of the usual anti-periodic boundary condition. This is what one would do at 
weak coupling. It was called the "B-formalism" in [16]. The other way is to keep 
anti-periodic boundary conditions, but instead use a shifted Hamiltonian, 



This coupling can be viewed as a coupling to a gauge field for which A = /i^. We 
will refer to this as the "A-formalism" (in which we will temporarily put tildes on all 
objects, as above, for clarity). For a non-anomalous symmetry, these two formalisms 
are equivalent, and one can go from the B-formalism to the A-formalism using a 
gauge transformation involving the external gauge field, with parameter 9a = —fiAt, 
relating the gauge fields in the two formalisms according to 




(3.27) 



H = H - haQl + HaQr ■ 




A/j, = Afj, + d„9 A ■ 



(3.29) 
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In terms of our holographic picture, this gauge transformation acts directly on 
the 5d gauge field, with a parameter Qa(z) that is ^-dependent, 



Qa{z) = tg A (z) 



g A (z -> ±00) = ±fi A . 



(3.30) 



which again acts on the gauge field as 4 



A 



■rn 



Am + <9 m G>A • 



(3.32) 



The difference between the two formalism can thus be formulated in a clear holo- 
graphic language as well: in the A-formalism the ansatz has Aq asymptoting to the 
chemical potential, whereas in the B-formalism the ansatz instead has this chemical 
potential stored in the A z component. The chemical potential is then best written 



which is nicely gauge invariant and independent of the formalism used. 

This is all clear and unambiguous when the symmetry is non- anomalous. How- 
ever, in the presence of an anomaly, one cannot pass from the formalism defined 
by (3.27) to that defined by (3.28). This is what happens in our system: as we 
discussed above, the 5d action after the anomaly correction (2.20) is not invariant 
under a chiral gauge transformation, 



The anomaly implies that one of the two formalisms is incorrect. 

The point of view of [16] is that there are strong indications that the B-formalism 
is the correct one in field theory. If one insists on computing with untwisted fermions, 
one needs to perform a gauge transformation, which not just introduces the chemical 
potential into A Q , but also modifies the action and Hamiltonian to correct for the 
fact that the action is not gauge invariant. To be precise, when one uses untwisted 
fermions, the action that one should use is the gauge-transformed action, which 
differs from the original one by the anomaly. This was called the "A'-formalism" 
(note the prime) in [16]. 

Let us see how this logic works for the Sakai-Sugimoto system under considera- 
tion here. The idea is thus that if we want to identify the chemical potential with 

4 Thc boundary gauge transformation parameter 9a{%) is related to the 5d parameter by 



as 




(3.33) 




:S[A e ] + S e [A e ]. (3.34) 



6a(x) — ^Qa(x, z — > ±00) 



MA*, 



(3.31) 



with the perhaps somewhat inconvenient signs following from our relation between the bulk and 
boundary gauge fields (2.9). 



the asymptotic value of Ao, we should be working not with S[Ae] but rather with 
S[Ai] + S e [Ai\. In order to compute the currents, we need to compute the variation 
of the 0a, which takes the form 

5S^ = ae m£npq J d i xdzd m \3(d n Q A )F^ q 5Aj + (d n Q A )^ q 5Af\ . (3.35) 

The contribution of the Q A term to the holographic currents can be obtained using 
the dictionary 



We then obtain 



dAfM [x, z = oo) 
AJ V = -6ae^{d v e A )F^ = Ga^e^F^ 

(3.37) 

AJ» = -2ae^{d v d A )F% = 2afi A e^F* . 

These expressions are independent of the particular function g A (z) which one chooses 
in (3.30). For our ansatz (3.6) with the boundary conditions (3.7) we get 

AJ° =0, 

AJ V = -12aBfi A , 



(3.38) 

v — L ^ u v ~ 11 



Aft = Aft = 



AJ° A = AJ\ = AJ\ = AJ\ = . 

This shows a very promising feature: any solution in this class will now exhibit the 
chiral magnetic effect, as there is a non- vanishing J v component. 

Unfortunately, we will see later that things are more subtle at the level of the 
Hamiltonian. There are two main problems when writing down the A'-formalism 
for things more complicated than the currents. Firstly, any bulk quantities such as 
the Hamiltonian will typically depend on g A (z), not just on its asymptotic values. 
Secondly, it will turn out that even for a 'natural' choice of g A (z), for instance 
g A (z) = —f A (z), the new Hamiltonian has the property that it does not lead to 
a minimum for non-homogeneous configurations. To see this requires some more 
details about this condensate, but let us here already present the expression for the 
corrected Hamiltonian. It can be written as 

H Tot {A e , B A ) = H Bulk (A e ) + H Bdy (A e ) + H^\A e , G A ) (3.39) 

where HBuik(Ai) and H-Qdy{At) are given by (2.27), (2.28) and the 0^ term is 

H^\A e ,Q A ) = aj d 3 xdz{e 0abcd (d b Q A ) [3 F^ d (d A^) + F*(d A*)} 

- \e im ^(d e e A )(3F^ q + J^)} (3.40) 
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For our ansatz (3.6) with the boundary conditions (3.7) the theta term in the Hamil- 
tonian take the form 



H^(Ai,Q A ) = ^a j d 3 xdz {d Q Q A ) 



-g^o / dzd z g A 



= -g^o J dzd z g A 
Specialising to g A (z) = —fA{z) we get 



3Bd z a v + -d z {3h 2 v + h 2 A ) 



k 



3Ba v + '-{3h 2 v + h 2 A ) 
3(&y + \h\) + \h\ 



(3.41) 



dz 



k 

-i 

2 



k 

bv + '^h 2 v + -h\ 



Kbv + \h 2 v ) + k -h\ 



-m pfi A . (3.42) 



In the homogeneous case we get 



H^\A E ,Q A ) 



-'HqC\ 



sinh(7ri?) 
B 



+ 7T 



^H Q pfi A 



(3.43) 



This boundary term can have drastic consequences for the phase structure of the 
theory: we will see in section 3.6 that it disfavours non-homogeneous configurations. 

We should emphasise that the procedure for introducing the A'-formalism in the 
holographic context is, unfortunately, rather ambiguous. This is essentially because 
the required boundary condition on the gauge transformation parameter Q A (z), given 
in (3.31), does not uniquely specify the behaviour of Q A (z) in the bulk. In contrast, 
this kind of ambiguity does not appear in field theory, as gauge transformation which 
untwists fermions, and "moves" chemical potential into temporal component of the 
gauge potential, is unique. 

Instead of using the A'-formalism one could consider doing the holographic com- 
putation directly in the B-formalism. This would require writing the solution in the 
A = gauge (the fermions are not directly accessible so it is unclear whether ad- 
ditional changes are required to implement the twisting). This does lead to a CME, 
but the Hamiltonian again turns out to have no minimum for non-homogeneous 
configurations. Some details are given in appendix A. 2. 



3.5 Perturbative stability analysis of the homogeneous solution 

In this section we will perturbatively analyse the stability of the homogeneous solu- 
tion (3.24), in order to show that this configuration is unstable and wants to "decay" 
to a non-homogeneous solution (3.6) (whose explicit form will be found later). Our 
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analysis here is a revision of the work done in [11], but our findings differ from theirs 
in an important way which is crucial for the remainder of our paper. 

Our starting point is given by the equations of motion linearised around the con- 
figuration (3.24). Following [11], we will look for fluctuations of the modes transverse 
to the direction of the external field, i.e. along (Ai,A 2 ), since these should lead to 
the formation of a chiral spiral. These fluctuations are given by 

5Ai = 5Ai{u, k)e- iut+lkxi (i = 2, 3) . (3.44) 

The equations of motion for the fluctuations (5A 1 ,5A 2 ) are coupled, but they diag- 
onalise in the complex basis 

<L4 (±) = 5A 2 ± i6A 3 , (3.45) 

where the equations become 

iC V V - k 2 )5A^ + Ml K d z (K z dM ±] ) ± ^-{kJF Qz + w7u)M (±) = . (3.46) 

Here T\ z , T§ z are field strengths of the background homogeneous solution (3.19), and 
we will express all results in terms of the constants Ca and Cb (instead of \xa and j) 
in order enable a simpler comparison with the results of [11]. 

Given values of Ca and Cb, we numerically solve equation (3.46). As usual 
in perturbation theory, solutions with real u represent fluctuations that are stable, 
while those for which u has a positive imaginary part correspond to instabilities, 
since they are exponentially growing in time. Fluctuations for which u = are 
marginal and have to be analysed in the full nonlinear theory in order to see if they 
correspond to unstable directions in configuration space. We will argue now that, 
while [1 1] has correctly identified the perturbatively unstable solutions with complex 
u, they have missed the marginally unstable modes, which are actually unstable in 
the full theory. In the following section we will then explicitly construct the new 
vacua corresponding to these marginal modes. 

Before presenting solutions to the equation (3.46), observe that this equation 
exhibits the symmetry t — > —t, u> — > —u, which means that all solutions will come 
in pairs (a;,— a;). Additionally, when solving this equation, one looks for normalis- 
able solutions, i.e. one looks for the solutions that behave as 5A^ ~ 1/z near the 
boundaries. 

Let us start by discussing the solutions with only Ca non-vanishing. At B = 
these correspond to solutions with j = 0. Samples of these solutions are shown in 
figure 1. We see that as \Ca\ is increased, the two branches of solutions come closer 
to the u = axis and then touch (the middle plot in figure 1). For all those values 
of \Ca\, given any momenta k, there is always real u solution. However, as \Ca\ is 
increased even more, the two branches of solutions separate along the /c-axis, so that 
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Q =-135542 C B =0. 



O =-2.02994, („=<) 



C A =-2.1031& C fl =0. 



Figure 1. Dispersion relation for small fluctuations at vanishing Cb and increasing values 
of Ca- In the last plot, the region between the two cu = modes corresponds to modes 
with positive imaginary oj, hence signalling proper instability. This follows Kim et al. [11]. 



C A =-1.95542, C B =-IM).")*xi) 



C^=-1.95542, C fl =-0.S' 



C,, =-1.95542 C s =-1.87312 





Figure 2. Dispersion relation at fixed value of Ca and increasing values of Cb- Modes 
with oj = occur for sufficiently large Cb, but in the region on the oj = axis between 
these two roots, there are no strictly unstable modes. Nevertheless, the oj could signal that 
there appears a new ground state. 

there is a region of momenta for which there are no real u> solutions (see the third 
plot of figure 1). For these "forbidden" values of the momenta, one can explicitly find 
solutions with complex u, which clearly signal a proper instability of the solution. 
These modes have previously been found in [11]. 

Let us now turn on \Cb\, while keeping |Ca| fixed. Samples of these solutions 
are shown in figure 2. We see that as |Cb| is increased, the two branches of solutions 
shift in the vertical direction, while the distance between them remains non-vanishing 
(which is the reason why this other branch is not visible in the second and third plots). 
For some value of \Cb\, the upper branch crosses the u = axis, and continues to 
go towards negative u. We thus see that for large enough \Cb\, marginal oj = 
modes are always present in the spectrum. In the next section, we will show that 
these modes, which were previously missed in the literature, are actually unstable 
once non-linearities are taken into account. 
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Figure 3. The physical picture (not to scale) that arises from the analysis of the dispersion 
relation. The plot displays the location of the physical modes at the extrema of the curves 
in figures 1 and 2, but now depicted in the complex u) plane. 



These findings are summarised in figure 3. Depicted there is the behaviour of 
the frequencies uj of the fluctuation modes at the 'extrema' of the dispersion relation 
plots, as a function of the two parameters Ca and Cg. For a generic value of Cb and 
Ca = 0, one root is always positive. As Ca is increased from zero, the positive root 
(red dot) moves towards the left, until it becomes zero. This defines the lower, solid 
curve of marginal modes in figure 3. Above this curve there is always a marginally 
stable mode in the spectrum, and thus potentially an instability. As Ca is increased 
further, both roots eventually develop a positive imaginary part, and we enter the 
region of strictly unstable modes (the blue shaded region). This area was previously 
discussed in [ll]. 5 



5 As we have already mentioned before, we disagree with their interpretation of Cb being related 
to the baryon chemical potential; sec also below. 
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3.6 The inhomogeneous solutions to the equations of motion 

Having established the location of the unstable and marginally stable modes in the 
fluctuation spectrum of the homogeneous solution, we would now like to explicitly 
construct the non-homogeneous solutions to which they are expected to decay, and 
study their properties. In our previous work [15] we have explicitly constructed a 
non-homogeneous solution in the presence of non-vanishing axial chemical potential, 
and in the absence of a magnetic field. We will therefore first discuss solutions with 
nonzero \ia and non-zero B, which are a natural generalisation of those considered 
before. We will then introduce a non-zero potential j as well and discuss how the 
solutions behave as a function of this parameter. 

The upshot of our previous analysis [15] was that for large enough axial chemical 
potential, larger than some critical value, a non-homogeneous solution is formed. 
The relation between the chemical potential and axial particle density was almost 
linear, similar to the homogeneous case. However, a particular density of particles 
was achieved for a smaller value of the chemical potential than in the homogeneous 
configuration. The wave number characterising the period of the non-homogeneous 
solution turned out to be very weakly sensitive to the actual value of the particle 
charge density in the canonical ensemble (or to the value of the chemical potential in 
the case of the grand canonical ensemble). The first things we would like to know is 
whether any of these characteristics change in the presence of an external magnetic 
field. 

Our starting point is the nonlinear system of equations (3.11), (3.12) and (3.13). 
We first observe that the function f(z) appears only in the first equation (3.11), and 
this equation can be directly integrated once we determine the functions b(z) and 
h(z), from the other two equations. 

Hence, we first need to solve equations (3.12) and (3.13). The parameters B and 
k in these equations correspond to the external magnetic field and the momentum of 
the transverse spiral, and are fixed at this stage. There are then four undetermined 
constants corresponding to non-normalisable and normalisable modes for each of the 
functions b and h. However, we solve equations (3.12), (3.13) requiring that the 
transverse spiral describes a normalisable mode, 

h{z) = ^ + • • • , (3.47) 

z 

i.e. we impose that the only external field is the magnetic field, and that there are 
no transverse external fields. 

In contrast, the function b(z) (or alternatively, a(z), see (3.14)), which describes 
the longitudinal field, does have a non-normalisable component. For example, on the 
positive side of the U-shaped brane it generically behaves as 

b{z) = b + - + --- (z^oo), (3.48) 
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Figure 4. Value of ho versus k at fixed B = 1.0, for increasing values of \Cb\- Increasing 
| (7b | makes the curve go up, indicating that a full new non- linear ground state exists for 
increasing \Cb\ and decreasing \Ca\- That is to say, the oj = fluctuation modes are the 
indicators of the transition to a new ground state, not the Im(u;) > modes. 



where the coefficients b and b\ are given by 

h i} ----- C.\ cosh ( -y - ) + Cb sinh ( ^ 



(3.49) 



6i = 5(7 B cosh f_J - sinh ^ 

However, just like for the homogeneous solution, this non-normalisable component 
of the solution corresponds to a gradient of the r( field and not to an external field 
(we will provide more evidence of this in section 3.6). Its unusual appearance as a 
non-normalisable mode is a consequence of our choice of the A z = gauge. 

We should note that at this stage we do not impose that b or h are of definite 
parity. We numerically solve the equation of motion using the shooting method, and 
to do this we only need to specify three undetermined constants h , Ca and Cb on 
one side of the U-shaped brane, as well as the parameter B and k in the equations. 
We solve equations for various values of constants and parameters, but keep only 
those for which h is normalisable. 

A first indication of the behaviour of the solutions can be obtained by looking at 
the effect of increasing Cb at non-zero value of B. Remember from section 3.5 that 
there exists a marginally stable mode in the spectrum, which follows the downward 
bending curve in figure 3. Along this curve there might be a decay to a new ground 
state. That this indeed happens is easily confirmed by looking for nonlinear solutions 
at non-zero Ca and Cb- A number of physical solutions is depicted in figure 4. They 
confirm that the marginal modes found in the previous section in fact correspond to 
true instabilities and new ground states. However, these figures say little about the 
actual physics, as neither Ca nor Cb are physical parameters. 




Figure 5. Data for solutions at Cb = 0, for the (arbitrary) value ll = —4.83. The second 
plot clearly shows that when B ^ 0, constant j solutions require a scan through the full 
parameter space {ho,k, Ca, Cb}- 

It is tempting to associate Cb with a baryon chemical potential, but as we have 
already mentioned, this is incorrect, as the parameter \i B completely decouples and 
does not influence the value of the Hamiltonian. Instead, the correct interpretation 
is that the physical parameters in our problem are fx a and j, which happen to be 
related in a non-trivial way to Ca and Cb- Families of solutions (parametrised by k) 
at fixed \xa and j lie on curves in the space spanned by Ca and C B , which only 
become straight lines at Cb = when B — 0. This is once more clearly visible in 
figure 5. This depicts a family of solutions at constant lla, parametrised by k, at 
Cb = 0. These clearly do not have a constant value of j. 

In order to find solutions with both lla and j fixed, we need to allow for a 
variation of both Ca and Cb as a function of k. Even with the magnetic field fixed 
to a particular value, this still means that the normalisable solutions lie on a curve in 
a four- dimensional parameter space spanned by {ho, k, Ca, Cb}- This makes a brute 
force scan computationally infeasible. Independent of the large dimensionality of this 
problem, we also found that the larger fx a cases require substantial computational 
time because of the fact that the asymptotic value h(— oo) varies rather strongly as 
a function of k and the other parameters. In other words, the valley of solutions is 
rather steep and the bottom of the valley, where h(z) is normalisable at z — > — oo, 
is difficult to trace. In this respect, it is useful to note that a solution written in 
C++ using odeint [29] and GSL [30] outperformed our Mathematica implementations 
by two to three orders of magnitude (!) in these computations. The details of the 
procedure which we followed can be found in appendix A. 3. 

We first focus on solutions at fixed and rather arbitrary value B = 1; the de- 
pendence on B will be discussed in the next section. We will mainly discuss j = 
solutions. Our numerical investigations of j ^ solutions show that all of these ac- 
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Figure 6. Left: Comparison of the energy of the non- homogeneous ground state (solid 
curve) with the energy of the homogeneous state (dashed curves) for j = 0, B = 1 and four 
values of \xa\ from top to bottom \xa = —5.5, —6.0, —6.5 and —7.5. Right: ho versus k, for 
the same values of the parameters (from bottom to top). 

tually have higher free energy, and are thus not real ground states of the system (see 
below). In figure 6 we display a set of configurations at constant [ia and vanishing j. 
Also depicted is the difference of the Hamiltonian of the homogeneous solution for 
this pair of fiA,j values, and the Hamiltonian of the non-homogeneous solution. 

As in earlier work without magnetic field [15, 27], there is a family of solutions 
parametrised by the wave momentum k, and the physical solution is the one for 
which the Hamiltonian is minimised. Perhaps somewhat surprisingly, for smaller 
values of the chemical potential (we have done computations up to about ha = —8) 
the momentum k g , s , at which the non-homogeneous ground state attains its minimum 
energy is only very weakly dependent on \ia- Our numerics for the range of ha up to 
—7.5 show that the ground state momentum for all these cases equals 0.83 to within 
less than one percent. Furthermore, this ground state momentum is also the same 
(to within numerical accuracy) as the ground state momentum for the B = case 
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Figure 7. The dependence of the Hamiltonian and the current J\ on the parameter j 
(showing j = 0, 0.1, 0.3 and 0.5 bottom to top on the left, top to bottom on the right). 
This shows that minimisation with respect to j drives the system to j = 0, both for the 
homogeneous and for the non-homogeneous state, and that this also drives the //-gradient 
sitting in j\ to zero. Under j — > —j the Hamiltonian is invariant while J\ changes sign. 



analysed in [15]. This is despite the fact that the actual solutions are quite different. 
It would be interesting to understand this behaviour better. 

If one evaluates the Hamiltonian of the A'-formalism (3.42), one finds that all 
non-homogeneous configurations always have higher energy than the corresponding 
homogeneous one at the same values of \ia and j. We take this as a strong sign 
that there is still something un-understood about the A'-formalism, as the analysis 
of section 3.5 clearly indicates a perturbative instability. It is in principle possible 
that we are not looking at the correct ansatz for the ground state, but we consider 
this unlikely. A similar statements holds for the B-formalism. 

The dependence of these solutions on j can also be computed, and is depicted 
in figure 7. From those plots one observes two things. Firstly, a minimisation of 
the Hamiltonian with respect to j drives the system to j = 0. Secondly, a non-zero 
j leads to a non- vanishing J\ current, which can be interpreted as an //-gradient. 
Together, these observations show that the preferred non-homogeneous state is one 
with a vanishing //-gradient. 
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Figure 8. The sign of h(z — > — oo), where red and orange regions denote minus and plus 
respectively, and grey regions correspond to solutions that diverge when z — > — oo. Nor- 
malisable non-homogeneous solutions sit at the boundary between red and orange regions. 
The branch that starts near k = 0.2 is the relevant one, as it has the lowest energy. It 
becomes increasingly difficult to find numerically as it is sandwiched in a steep valley as 
k increases. All plots are at Ca = —2.4, with B = 1.9, 2.1 and 2.3 respectively. Note the 
different axis scales on the first plot. 




Figure 9. Left: The increase of B from to 1.0 to 1.5 (bottom black to top orange 
curves) at Cb = and fixed Ca = —2.4 has the effect of raising the curves of normalisable 
solutions. This implies that a larger magnetic field requires a smaller value of \Ca\ to 
trigger an instability. Right: However, one cannot conclude anything about the effect at 
constant [xa from this, since this physical parameter varies over the curves in the figure on 
the left. 
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Solutions at fi A = -6.0, j = for D = 1.00, 1.20, 1.40, 1.59 (top to bottom) 



0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 




2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 



Figure 10. Solutions at constant fi = —5.0 (left) and \x = —6.0 (right) and j = 0, for 
various values of the magnetic field. This shows that a magnetic field suppresses the insta- 
bility to a non-homogeneous ground state, and there is a critical magnetic field B CT i t (l^A) 
above which the non- homogeneous solution ceases to exist. 



Figure 11. Results from a linearised analysis close to the critical magnetic field. Left: 
critical chemical potential (i.e. the potential for which a non-homogeneous ground state 
first develops) versus critical magnetic field (i.e. the field for which the solution disap- 
pears). Right: as determined earlier for the non-linear case, the parameter Ca goes to zero 
(asymptotically) as B increases. 

3.7 Critical magnetic field 

For larger magnitudes of the external magnetic field, the numerics becomes increas- 
ingly expensive as the valleys of solutions to the equations of motion rapidly become 
steeper and are more and more closely approached by regions in parameter space 
in which no regular solutions exist (i.e. regions in which h(z — > oo) diverges). See 
figure 8 for an impression. 

A first analysis which one can make is to simply scan for normalisable solutions 
at fixed Ca and Cb, for increasing values of B. This leads to plots like in figure 9. 
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Figure 12. Dependence of the ground state momentum fc g . s . of the critical solution on 
the magnetic field or the chemical potential, computed in the linear approximation. The 
non- linear solutions which we computed cover the regime up to about B c = 2. One should 
keep in mind that the linear approximation will break down for sufficiently large values of 
B c . 

One observes that a larger magnetic field requires a smaller \Ca\ for the chiral spiral 
solution to exist. In [11] this effect was seen at the level of the perturbative insta- 
bility analysis, and led the authors to conclude that a magnetic field enhances the 
instability. However, this conclusion is premature and actually incorrect. What one 
needs to do is to analyse the effect of the magnetic field on solutions at fixed value 
of /I a and j, not at fixed values of the unphysical parameters Ca and (7b . 

If one does this more elaborate analysis, the conclusion is actually opposite: the 
magnetic field tends to stabilise the homogeneous solution, and there exists a crit- 
ical B crit (fiA) above which the non-homogeneous solution ceases to exist. This can 
be seen in figure 10, where solutions at fixed fx a and j are depicted. By increas- 
ing B sufficiently slowly while keeping the other physical parameters fixed, we can 
determine this -B cr i t in numerical form. 

It is difficult to get a good picture of B CTit versus fiA, as the numerics become 
expensive, for reasons we have mentioned earlier. However, we can make use of 
the fact that near the critical magnetic field, the parameter ho is small. Assuming 
that this implies that the entire function h(z) is small, we can then use a linear 
approximation to the equations of motion, which is much easier to solve. The result 
of this linear analysis is depicted in figures 11 and 12. The latter shows that the 
ground state momentum k g s is actually not as flat as the non-linear computations 
suggests. 

For larger B, one should keep in mind that the results may be invalidated for a 
variety of reasons. Firstly, the linear approximation may break down, as smallness of 
ho does not necessarily imply smallness of the full function h(z). Secondly, we have 
seen that the valley near the non-linear solutions becomes very steep for large values 
of B, and hence the linear solution may deviate quite strongly from the non-linear 



one for large B. Thirdly, when B is large, DBI corrections to the equations of motion 
may become relevant, as higher powers of the field strength are no longer necessarily 
small. For these reasons, one should be careful interpreting the large B c regime of 
figure 12. 

4 Conclusions and open questions 

We have analysed in detail the instability of the Sakai-Sugimoto model in the presence 
of a chiral chemical potential [ia, an //-gradient j and a magnetic field B. We have 
shown that the presence of marginally stable modes (overlooked in [11]) is a signal 
for decay towards a non-homogeneous chiral spiral like ground state. Minimising 
the Hamiltonian on a constant fiA, constant j curve leads to a unique ground state 
momentum k g , s , (see figure 6), which for small B is only weakly dependent on fiA- 
Increasing the magnetic field to sufficiently large values suppresses the instability 
and drives the system back to the homogeneous phase (figure 12). This result may 
be related to the effective reduction of the Landau levels in a strong magnetic field 
(see [31]). 

A linear analysis suggests that the ground state momentum may have a non- 
trivial dependence on the chemical potential, and might be compatible with a linear 
scaling at large ^a (as in the Deryagin-Grigoriev-Rubakov non-homogeneous large- N c 
QCD ground state [32]). However, substantial additional work is necessary to deter- 
mine whether this is indeed happening in the full non-linear theory. 

We have found that a recent proposal for the correction of the currents [13, 16], 
while capable of producing the chiral magnetic effect both in the homogeneous and 
non-homogeneous ground states, is incomplete. It is a) not unique and b) leads to a 
Hamiltonian which does not prefer the non-homogeneous ground state. We emphasise 
that this is a problem with the currents and boundary terms in the Hamiltonian, as 
the perturbative instability analysis is not affected by these corrections, and neither 
are the non-linear solutions. 
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A Appendix: Technical details 



A.l The Hamiltonian after the anomaly corrections 



Here we will give a detailed derivation of the Hamiltonian for the Sakai-Sugimoto 
system taking special care about surface terms that are present, as well as taking 
into account the corrections to the action coming from the anomaly. 

After anomaly in the system have been corrected by inclusion of extra terms in 
the action, the full Langrangian density can be written as 



C = — m 



-9 



pOa pV , pab pV , -pOa pA , pab pA 
V Oct * 2 V- r a6" l "* r A* r 0a"T ^ A ab 
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Oabcd 



A^I b ^ + ^ b ^)+^Tl^ + T A a ^ . (A.l) 
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The conjugate momenta associated with the vector and axial gauge fields take 
the form 

f)T 



ll A 



d{d Al) 
dl 
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cd ■ 



d(d Ai) 

Hence, the on-shell Hamiltonian takes the form H = H Bu i k + H Bdy 

1 



#Bulk — K J 
H^dy = j 



d xdzy— <7 



pOa p I pab p 
~J~ J-Qa + ~J~ J~ ab 



d 3 :rd,2 d„ 



W V A V Q + W A A A 



(A.3) 
(A.4) 



where we have used the gauge field equations for the time component of gauge po- 
tential (generalised Gauss law). 

For our inhomogeneous ansatz the conjugate momenta simplify to 



fl z v = 2n^g zz g m d z f v - Qa(a A B + kh A h v ) = , 

ITy = —6a [aA(d z hv) — h A (d z ay)} [sin(/ca;i)a;2 + cos (/err 1)2:3] , 

fl z A = 2Ky/^g~g zz g 00 d z f A - 2akh\ = 6aBa v + 3akh 2 v + akh\ - 2np , 

IU = —2a [aA(d z hji) — h A (d z aA)] [sin(A;a;i)a; 2 + cos(kxi)x 3 ] . 



(A.5) 



Using this result and the boundary conditions we have imposed on the fields, one 
gets the boundary term of the Hamiltonian 



#Bd y = -2/cp J d 3 xdzd z f A = 4kV x pii a = -V x J A [i A . 



(A.6) 
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The bulk term simplifies for our ansatz to 



Bulk 



Kid xdz 



99 zz g 00 (d z ff 



+ v <n<rr [b 2 + k^} } . 



-ggPf* [{d z hf + (d z a) 2 ] (A.7) 

(A.8) 



In terms of the variables (3.16) the bulk and boundary terms take the form 

2" 
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+ K Z 
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(A.9) 

(A.10) 
(A.ll) 



where H Q = M^ K V x \ 3 N c /(8ir 2 ), b = Ba - p and H Div = n B 2 Jdz K z 1/Z . For the 
homogeneous case we get simple analytic result 
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^Buik = 2H B coth(|B)(/& + j 2 ) + H Di - 

n (dl + c 2 B ) shlh ^ 
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Div 



H Bdy = -m Bcoth(^B)fa = -2H Ci Sinh(57r) 



B 



(A.12) 
(A.13) 



A. 2 The B-formalism ansatz 



In formalism B the time component of the chiral gauge field is zero at the boundary. 
An ansatz that preserves the field strengths (and hence the equations of motion) is 
given by 



A 1 



-t d-J\ 



At 



-t d z f A 



At = 



B 



A v = — [-x 3 x 2 + x 2 x 3 ] + h v (z) [cos(kxi)x 2 - sin(/cxi)x 3 ] + a v (z)x 1 

•Aa = (Ia(z)xi + Iia(z) [cos(kxi)X2 — sin(A;xi)x3] , 
which should be used together with the boundary conditions 

A v z {z -> ±oo) = Af(z -> ±oo) = , 
hy{z — > ±oo) = Ha(z — > ±oo) = , 
ay{z — »■ ±oo) = a^(z — > ±oo) = =Fj . 



(A. 14) 



(A.15) 
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The components of the vector and axial currents take the same form as those in the 
A-formalism, given in (3.25), with the exception of J v which reads 

J v = -Ak lim \y/ = gg" x g xa d z a v ] = \2aB lim f A = -12aBfx A . (A. 16) 

This in particular implies that there is a non-zero chiral magnetic effect in formal- 
ism B. 

A. 3 Scanning parameter space for solutions 

We recall from the main text that we aim for the minimisation of the Hamiltonian 
as a function of the condensate momentum k, on curves which have fixed values of 
ix a and j. Unfortunately, the latter two parameters are determined only indirectly, 
after a solution has been found; their dependence on the shooting parameters h , Ca 
and Cb is not known analytically. Minimising the Hamiltonian over curves at fixed 
value of Ca and Cb would in general not be equivalent (i.e. wrong). 

One could in principle make a fine-grained lattice scan for solutions in the pa- 
rameter space spanned by {ho, k, Ca, Cb}, and then consider only those points which 
have a certain fixed value of \x A and j. This, however, is computationally extremely 
expensive, and wastes a lot of time on regions of parameter space which will never 
be used. We therefore follow a different approach. 

The general problem which we need to solve for an efficient determination of 
the required curves is the following. We have a <i-dimensional parameter space of 
configurations, and physical solutions are those in which d — 1 functions of the con- 
figuration are vanishing simultaneously. One of these functions is always the value 
of the function h(z) evaluated at — oo, as this imposes normalisability. Among the 
other functions we have, for instance, the value of fx a for this configuration minus 
some fixed reference value, in case we want to scan for solutions at this fixed value of 
the chemical potential. This 'isocurve tracing' problem is most easily implemented 
by starting from one known point on which the functions are simultaneously zero, 
and then using a version of the Newton-Raphson method to find the location of a 
neighbouring simultaneous zero. 6 If we have n lattice points in every direction of 
parameter space, this brings the computational cost down from being order 0(n d ) 
to the much more tractable 0(ri). 

We start by scanning for physical solutions at Cb = and some fixed value of 
Ca, which is a two-dimensional problem in {ho, k} space and yields curves like those 
in figure 5. On these curves neither \x A nor j will in general be constant. However, 
one can choose one point on such a solution curve, and use this as a seed point p seed 
for the isocurve tracing described above. Denoting the value of \x A and j by //^ ed 
and j' seed respectively, we then trace the common zero of the three functions h(—oo), 

6 General purpose Mathematica and C++ implementations of this isocurve tracing algorithm are 
available upon request from the authors. 
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1.0 



Figure 13. Typical curve of solutions, here at constant ha = —4.83 and constant j = 0.64 
(thick curve) and for B = 1. The orange/yellow surfaces show the projections onto the 
{A;, Ca} and {k,Cs} planes respectively. Only for vanishing magnetic does this curve lie 
completely inside the {k, Ca} plane. 

\i a — \± s< a A an d J ' ~ j seed m the four-dimensional parameter space {ho, k, Ca, Cb} (only 
when B = do these curves lie completely in the {ho,k,CA} subspace). A typical 
solution is depicted in figure 13. The only remaining problem with these curves is 
now that in general j seed ^ (the only such points on the curves in figure 5 are 
located on the endpoints of those curves, where the coefficient h goes to zero). 

In order to find curves at j = 0, we start again from p seed ; but first do an isocurve 
trace at fixed k in the {ho, Ca, Cb} space until we find a point at which j = 0. An 
example is given in figure 14. This point is then used as our new seed point p seed for 
the four-dimensional scan. Repeating the whole process using a seed point obtained 
for different initial values of Ca then produces a set of curves at j = for various 
constant values of \ia- 

On each of these curves we can now finally compute the Hamiltonian as a function 
of k, and find the value k m i n for which it is minimised. For a sample set of values of 
ji a the result is depicted in figure 6. 
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Figure 14. Typical isocurve trace, here at the sample value of ha = —4.83, starting from 
the seed point p seed with Cb = to the new seed point p seed (indicated by the dot) at 
which j = and both Ca and Cb are non-zero. 
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